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INTRODUCTION 


Electromechanical  transduction  taking  place  via  the  piezoelectric 
effect  is  characterized  phenomenologically  by  constitutive  equa¬ 
tions  that  relate  the  elastic  and  electric  variables.  These  equa¬ 
tions  take  a  variety  of  forms,  depending  upon  the  choice  of  inde¬ 
pendent  and  dependent  variables;  the  choice  is  normally  dictated  by 
the  application.  For  example,  piezoelectric  resonators  in  the  form 
of  thickness  mode  plates  are  most  easily  treated  using  the  isagric 
elastic  stiffnesses  [cE] ,  the  piezoelectric  stress  constants  [e] , 
and  the  dielectric  permittivities  at  constant  strain  [(ens)Sl. 

Various  measurement  techniques  yield  values  for  the  elements  of  a 
particular  coefficient  set  more  directly  than  those  of  another. 
The  coefficients  appearing  in  the  different  equation  sets  are, 
however,  interrelated,  so  that  once  any  one  complete  set  is  avail¬ 
able,  all  the  other  sets  of  elements  may  be  found.  The  most  accu¬ 
rate  and  precise  experimental  results  to  date  have  been  from  plate 
resonator  (resonance)  and  pulse-echo  (transit-time)  measurements. 
From  the  [cE],  [e],  and  [(eps)S]  matrices  determined  therefrom, 
those  matrices  representing  material  properties  expressed  in  the 
other  alternative  forms  may  be  calculated. 

Electrooptical  applications  are  becoming  increasingly  important. 
So  also  are  treatments  of  piezoelectric  and  ferroelectric  phenomena 
from  the  standpoint  of  molecular  interactions.  In  both  of  these 
cases  the  constitutive  equations  using  polarization  as  the  inde¬ 
pendent  electrical  variable,  rather  than  either  electric  intensity 
or  displacement,  assume  greater  importance  than  the  sets  tradition¬ 
ally  used  for  transducer,  signal  processing,  and  resonator  applica¬ 
tions  . 

In  this  report  we  give  the  complete  sets  of  linear  constitutive 
equations  relating  elastic  and  electric  fields.  For  each  equation 
set  the  numerical  values  are  computed  for  lithium  tantalate,  a 
highly  piezoelectric  ferroelectric  with  temperature-compensated 
orientations,  from  the  measured  [cE],  [e],  and  [(eps)S]  values  of 
Smith  and  Welsh  (Ref.  1).  Coupling  to  the  thermal  field  is  ne¬ 
glected.  Rationalized  mks  units  are  used  throughout. 


CONSTITUTIVE  EQUATION  SETS 

Symbols  and  units  for  the  quantities  employed  are  given  in  Table  1. 
In  terms  of  these,  six  constitutive  equation  sets  are  used.  Of 
these,  electric  intensity,  dielectric  displacement,  and  polariza¬ 
tion  each  appear  in  two  sets  as  an  independent  variable.  The  sets 
are,  in  compressed  matrix  notation,  as  follows.  A  prime  denotes 
transpose;  [I]  is  the  unit  matrix. 


I.  The  Piezoelectric  Stress  Constant  Set 

[T]  =  [cE]  [S]  -  [e]'  [E]  (  1) 

[D]  =  [  e]  [S]  +  [ (eps)Sj  [E]  (  2) 


TABLE  1.  SYMBOLS,  UNITS,  AND  DEFINITIONS. 


QUANTITY 

UNIT 

SYMBOL/DEFINITION 

Elastic  stress 

Elastic  strain 

N/m2 

[T] 

[S] 

Electric  intensity 

V/m 

[E] 

Dielectric  displacement 

C/m2 

[D] 

Dielectric  polarization 

C/m2 

[P] 

Elastic  compliance  at 
constant  [E],  [D] ,  [P] 

m2/N 

[CE],  [CD],  [CP] 

Elastic  stiffness  at 
constant  [E]  ,  [D] ,  [P] 

N/m2 

[sE],  [sD],  [sP] 

Dielectric  permittivity 
at  constant  [T] ,  [S] 

Dielectric  constant, 
relative,  at 
constant  [T] ,  [S] 

F/m 

[ (eps) T] ,  [ (eps) S] 

[(  Kr) T] ,  [(  Kr) S] 

=[ ( eps ) T ] / ( eps ) o , 
t (eps) S]/ (eps) o 

Dielectric  impermeability 
at  constant  [T] ,  [S] 

m/F 

[(bet)T]  [fbet)S] 

- [ (eps) T]  "I  , 

[ (eps) S ]  ( 

Dielectric  impermeability, 
relative,  at  constant 
[T],  [S] 

[ (betr) T] ,  [ (betr) S ] 

=  [ (bet) T] * (eps) o, 

[ (bet) S ] * (eps) o 
=[(  Kr) T]  "I  , 

[(  Kr) S ] 

Dielectric  susceptibility 
at  constant  [T] ,  [S] 

F/m 

[ (chi) T] ,  [ (chi) S] 
=[(  Kr) T-I ] * (eps) o, 

[ (  Kr) S-I ] * (eps) o 

Dielectric  susceptibility, 
relative,  at  constant 
[T],  [S] 

[ (chir) T] ,  [ (chir) S] 
=[ (chi) T]/ (eps) o, 

[ (chi) S]/ (eps) o 

Reciprocal  dielectric 
susceptibility  at 
constant  [T] ,  [S ] 

m/F 

[ (zet)T],  [fzet)S] 

= [ ( chi) T]  , 

[ ( chi ) S ]  (-1) 

Reciprocal  dielectric 
susceptibility,  relative, 
at  constant  [T] ,  [S] 

[ ( zetr) T] ,  [ ( zetr) S ] 
=[ ( zet ) T] * (eps) o, 

[  ( zet) S ] * (eps) o 

Piezoelectric  stress 

C/m2 

[e] 

Piezoelectric  stress 
constant 
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TABLE  1.  SYMBOLS,  UNITS,  AND  DEFINITIONS,  (continued) 


QUANTITY 

UNIT 

SYMBOL/DEFINITION 

Piezoelectric 

coefficient 

strain 

m/V  =  C/N 

[d] 

Piezoelectric 

modulus 

stress 

N/C  =  V/m 

[h] 

Piezoelectric 

constant 

strain 

m2/C 

[g] 

Piezoelectric 

modulus 

polai ization 

V/m  =  N/C 

[a] 

Piezoelectric 

constant 

polarization 

m2/C 

[b] 

Note:  Square  brackets,  sic :  [  ],  denote  matrices. 
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II.  The  Piezoelectric  Strain  Coefficient  Set 

[SJ  =  [sE]  [T]  +  [d]'  [E]  (  3) 

[D]  =  [  d]  [T]  +  [ (eps)T]  [E]  (  4) 

III.  The  Piezoelectric  Stress  Modulus  Set 

[T]  =  [CD]  [S]  -  £h] •  [D]  (  5) 

[E]  =  -[  h]  [S]  +  [  (bet)  S ]  [D]  (  6) 

IV.  The  Piezoelectric  Strain  Constant  Set 

[S]  =  [SD]  [T]  +  [g]'  [D]  (  7) 

[E]  -  -[  g]  [T]  +  [ (bet) T]  [D]  (  8) 

V.  The  Piezoelectric  Polarization  Modulus  Set 

[T]  =  [cP]  [S]  -  [a]'  [P]  (  9) 

[E]  =  -[  a]  [S]  4-  [  ( zet)  S  ]  [P]  (10) 

VI.  The  Piezoelectric  Polarization  Constant  Set 

[S]  =  [SP]  [T]  +  [b]'  [P]  (11) 

[E]  -  -[  L]  [T]  +  [izet)T]  [P]  (12) 


The  electric  variables  are  connected  by  the  relation 

[D]  =  (eps)o  *  [E]  +  [P]  (13) 

where  (eps)o  is  the  permittivity  of  free  space,  defined  by 

(eps)o  *  (mu)o  *  (c)  *  (c)  =  1  ;  (14) 

(mu)o  is  the  permeability  of  free  space,  equal,  by  definition,  to  4 
PI  *  10(-7),  and  (c)  is  the  velocity  of  light  in  vacuo  and,  also  by 
definition,  is  equal  exactly  to  2.99792458  x  108  m/s. 


From  (13)  the  expressions  for  the  remaining  electric  variables 
associated,  respectively,  with  the  six  equation  sets  (  1)  to  (12) 
may  be  found: 


[P]  = 

[  e]  [S] 

+  [ (chi) S ] 

[E] 

(15) 

[P]  = 

(  d]  [T] 

+  [ (chi) T] 

[E] 

(16) 

[P]  = 

(eps)o  * 

[  h]  [S]  + 

[I  - 

(eps)o  * 

(bet) S] 

[  D 1 

(17) 

[P]  = 

(eps)o  * 

[  g]  [T]  + 

[I  " 

(eps)o  * 

(bet) T] 

[D] 

(18) 
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[D] 

=  - ( eps ) o  * 

[  a] 

[S]  + 

[I  + 

(eps)o  * 

( zet) S ] 

[P] 

(19) 

[D] 

=  - (eps) o  * 

[  b] 

[T]  + 

[I  + 

(eps)o  * 

( zet) T ] 

[P] 

(20) 

RELATIONS  AMONG  MATERIAL  CONSTANTS 

The  material  constants  are  interrelated  by  the  following  formulas: 

[cX]  [sX]  =  [ ( eps ) Y ]  [ (bet) Y ]  =  [I]  (21) 

[ (chi) Y]  [ ( zet ) Y ]  =  [(  Kr ) Y  -  (chir)Y]  =  [I]  (22) 

In  (21)  and  (22),  X  =  E,  D,  or  P  and  Y  =  T  or  S . 

(CD]  -  [cE]  =  [h] ’  [e]  =  [e] '  [h] 

=  [h] •  [ ( eps) S ]  [h]  -  [e] *  [ (bet) S ]  [e] 

=  [a]'  [e  -  h  *  (eps)o]  =  [e  -  h  *  (eps)o]'  [a]  (23) 

[cP]  -  [cD]  =  [h]'  [a]  *  (eps)o  =  [a]'  [h]  *  (eps)o 

=  [h]'  [(eps)S]  [ ( zet) S ]  [h]  *  (eps)o 

=  [a]'  [ (bet) S ]  [ (chi) S ]  [a]  *  (eps)o 

=  [a  -  h] '  [e]  =  [e] •  [a  -  h]  (24) 

[CP]  -  [cE]  -  [a] '  [e]  =  [e] '  [a] 

=  [a]'  [ (chi ) S ]  [a]  =  [e] •  ((zet)S]  [e] 

=  [h] '  [e  +  a  *  (eps)o]  =  ie  +  a  *  (epc)o]'  [h]  (25) 

[sE]  -  [sD]  =  [d] 1  [g]  =  [g] '  [d] 

=  [d] •  [ (bet) t]  [d]  =  [g] •  [ ( eps) t ]  [g] 

=  [ b ] •  [d  -  g  *  (eps) o]  =  [d  -  g  *  (eps)o]'  [b]  (26) 

[sD]  -  [sP]  =  [b]'  [g]  *  (eps)o  =  [g] ’  [b]  *  (eps)o 

=  [g] '  [ (eps) T]  [ ( zet) T ]  [g]  *  (eps)o 

=  [ b ] '  [ ( bet ) T ]  [ (chi ) T]  [b]  *  (eps)o 

=  [b  -  g] '  [d]  =  [d] *  [b  -  g]  (27) 

[sE]  -  [sP]  =  [b] '  [d]  =  [d] '  [b] 

=  [b] •  ( ( ch i ) T ]  [b]  =  [d]*  [ ( zet ) T ]  [d] 

=  [g] '  [d  +  b  *  (eps)o]  =  [d  +  b  *  (eps)o] '  [g]  (28) 
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[ ( zet) S ]  -  [ ( zet) T]  =  [b]  [a]'  =  [a]  [b] ' 

=  [b]  [cP]  [b] '  =  [a]  [sP]  [a] •  (29) 

[ (chi) T]  -  [ (chi) S ]  =  [(eps)T]  -  [(eps)S] 

=  [e]  [d]'  =  [d]  [e] • 

=  [d]  [cE]  [d] •  =  [e]  [sE]  [e] •  (30) 

[ ( bet ) S ]  -  [ (bet) T]  =  [h]  [g]'  =  [g]  [h] • 

=  [g]  [CD]  eg] '  =  [h]  [SD]  [h] •  (31) 

[e]  =  [d]  [CE]  =  [ ( eps ) S ]  [h]  =  [(chi)S)  [a]  (32) 

[d]  =  [e]  [sE]  =  [ (eps) T]  [g]  =  [(chi)T]  [b]  (33) 

[h]  =  [g]  [CD]  =  [ (bet) S]  [e]  =  [(chi)S]  [(bet)S]  [a] 

=  [I  -  (bet) S  *  (eps) o]  [a]  (34) 

[g]  =  [h]  [sD]  =  [ (bet) T]  [d]  =  [(chi)T]  [(bet)T]  [b] 

=  [I  -  (bet) T  *  (eps) o]  [b]  (35) 

[a]  =  [b]  [cP]  =  [ ( zet) S ]  [e]  =  [(eps)S]  [(zet)S]  [h] 

=  [I  +  (zet)S  *  (eps)o]  [h]  (36) 

[b]  =  [a]  [sP]  =  [ ( zet ) T ]  [d]  =  [ (eps) T]  [(zet)T]  [g] 

=  [I  +  ( zet) T  *  (eps) o]  [g]  (37) 

Some  alternative  relations  are  the  following: 

[a  -  h]  =  [(zet)S]  [h]  *  (eps)o 

=  [(bet)S]  [a]  *  (eps)o  (38) 

[b  -  g]  =  [ ( zet) T ]  [g]  *  (eps)o 

=  [ ( bet ) T ]  [b]  *  (eps)o  (39) 

[e  +  a  *  (eps)o]  =-  [(eps)S]  [a]  (40) 

[d  +  b  *  (eps)o]  =  [(eps)T]  [b]  (41) 

[e  -  h  *  (eps)o]  =  [(chi)S]  [h]  (42) 

[d  -  g  *  (eps) o]  =  [(chi)T]  [g]  (43) 


Equations  (21)  to  (43)  result  from  equating  like  dependent  varia¬ 
bles  in  pairs  selected  from  equations  (  1)  to  (12)  and  (15)  to 


pair  yields  one  equation  in  three  variables,  one  mechanical  and  two 
electrical,  or  vice  versa.  Two  other  equations  exist,  again  from 
(1)  to  (12)  and  (15)  to  (20),  that  contain  the  same  three  variables 
found  in  each  paired  equation.  One  of  these  auxiliary  equations  is 
used  to  eliminate  one  of  the  two  variables  of  the  same  kind;  the 
result  is  one  equation  in  two  variables,  one  electrical  and  one 
mechanical.  These  are  now  independent  variables,  so  the  coeffi¬ 
cients  must  vanish;  two  relations  between  the  material  coefficients 
result.  As  an  example,  (  3)  and  (  7)  both  have  [S]  as  dependent 
variable.  Equating  them  produces  one  relation  in  [T]  ,  [E],  and 
[D);  one  of  the  electrical  variables  must  be  eliminated.  This  is 
done  by  using  either  (  4)  or  (  8);  each  contains  the  same  three 
variables.  If  (  8)  is  used  to  eliminate  [E],  one  obtains  [sE  -  d' 
g  -  sD]  [T]  =  [d'  (bet) T  -  g ' ]  [D].  Therefore,  [sE]  -  [sD]  =  [d] ' 
[g]  and  [g]  =  [(bet)T]  [dj.  Use  of  (  4)  instead  of  (  8)  leads  to 
the  equations  [sE]  -  [sD]  =  [g] '  [d]  and  [d]  =  [(eps)T]  [g] .  There 
are  36  pairs,  six  each  equating  [S]  and  [T] ,  and  eight  each  equat¬ 
ing  [E],  [D],  and  [P].  The  72  relations  contain  many  redundancies. 
Relations  between  the  elastic,  piezoelectric,  and  dielectric  con¬ 
stants  are  shown  schematically  in  Tables  2  and  3. 


CALCULATION  SEQUENCE 


Using  as  input  [cF],  [e],  and  [(eps)S],  one  may  compute  the  remain¬ 
ing  quantities  in  a  variety  of  ways.  The  following  sequence  is 
typical : 


[sE]  =  [cE] 

(44) 

[ ( bet ) S ]  =  [ (eps)S] 

(-1) 

(45) 

[d]  =  [e]  [sE] 

(46) 

[h]  =  [ (bet) S ]  [e] 

(47) 

[ (eps) T]  -  [ (eps) S ] 

=  [e]  [d )  ' 

(48) 

[ ( eps ) T ]  =  [ (eps)S] 

+  [e]  [d]1 

(49) 

[ (bet) T]  =  f (eps) T] 

(-D 

(50) 

[cD]  -  [cE]  =  [ e • ' 

[h] 

(51) 

[cD]  -  [CE]  4  [e] ' 

[h] 

(52) 

[g]  =  [ ( bet ) T ]  [d] 

(53) 

[sE]  -  [sD]  =  [d] ' 

rU] 

(54) 

[sD]  -  [sE]  -  [d] ' 

[q] 

(55) 

[  (betr) S]  =  [ (bet) S 

]  *  (eps)o 

(56) 

[  (  zetr ) S ]  =  [  ( betr) 

S]  (I  -  ( bet  r ) S ]  <  ~  1 ) 

(57) 

(58) 


[  ( zet ) S ]  =  [(zetr)i  /  (eps)o 
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TABLE  2.  RELATIONS  AMONG  MATERIAL  CONSTANTS. 
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FURTHER  RELATIONS  AMONG  MATERIAL  CONSTANTS. 


[ (betr) T]  =  [ (bet) T]  *  (eps)o 

(59) 

[(zetr)T]  =  [ (betrT]  [I  -  (betr)T] 

(-D 

(60) 

( (zet)T  =  [(zetr)T] 

/  (eps)o 

(61) 

[ (chi) S]  =  [ ( zet) S ] 

(-1) 

(62) 

[ (chi) T]  =  [ (zet) T] 

(-1) 

(63) 

[a]  =  [ ( zet) S ]  [e] 

(64) 

[b]  =  [ ( zet) T]  [d] 

(65) 

[cP]  -  [cE]  =  [e] • 

[a] 

(66) 

(cP]  =  [cE]  +  [e]' 

[a] 

(67) 

[cP]  -  [CD]  =  [a] ' 

[h]  *  (eps)o 

(68) 

[sE]  -  [sP]  =  [d] ' 

[b] 

(69) 

[sP]  =  [sE]  -  [d]' 

[b] 

(70) 

[sD]  -  [sP]  =  [g]' 

[b]  *  (eps)o 

(71) 

[ (bet) S ]  -  [ (bet) T ] 

=  [h]  [g]' 

(72) 

[ ( ch i ) T ]  -  [ (chi) S ] 

=  [ (eps) T]  -  [ (eps) S ] 

(73) 

[ ( zet) S ]  -  [ ( zet) T] 

=  [a]  [b] • 

(74) 

A  number  of  these 

relations  are 

used  as  checks. 

For  example 

[(bet)S]  and  [(bet)T]  are  known  from  (45)  and  (50),  but  the  differ¬ 
ence  is  recomputed  in  (72) . 

EXPLICIT  FORMULAS  FOR  POINT  GROUP  3m 
Elastic : 

The  6x6  elastic  constant  portions  of  Tables  4  and  5  partition  into 
4x4  and  2x2  submatrices.  The  4x4  elastic  stiffness  and  compliance 
submatrices  are  interrelated  by  formulas  (75)  to  (93),  taken  from 
Cady  (Ref .  2 )  : 


A  = 

s33 

*  (sll  +  sl2)  - 

2 

* 

sl3 

* 

sl3 

(75) 

B  = 

s44 

*  (sll  -  sl2 )  - 

2 

★ 

sl4 

★ 

sl4 

(76) 

2  * 

ell 

=  s3 3  /  A  +  s4 4 

/ 

B 

(77) 

2  * 

cl2 

=  S3 3  /  A  -  s44 

/ 

B 

(78) 

Cl3 

=  - 

sl3  /  A  ;  cl4 

= 

- 

sl4 

/ 

B 

(79a),  (79b) 

(79) 
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TABLE  4.  ELASTOPIEZODI ELECTRIC  MATRICES  FOR  POINT  GROUP  3m: 


THE 

[e]  / 

[h] , 

AND  [a] 

SETS. 

11 

12 

13 

14 

00 

00 

]  00 

-22 

31 

] 

CE 

]  e' 

12 

11 

13 

-14 

00 

00 

)  00 

22 

31 

— 

-] - 

] 

e 

] (eps)S 

13 

13 

33 

00 

00 

00 

]  00 

00 

33 

14 

-14 

00 

44 

00 

00 

J 

]  00 

15 

00 

] 

CD 

]  h' 

00 

00 

00 

00 

44 

14 

]  15 

00 

00 

— 

■] - 

3 

h 

] (bet) S 

00 

00 

00 

00 

14 

66 

]  -22 

00 

00 

00 

00 

00 

00 

15 

-22 

J 

]  11 

00 

00 

] 

CP 

3  a1 

-22 

22 

00 

15 

00 

00 

]  00 

11 

00 

— 

■] - 

] 

a 

] ( zet) S 

31 

31 

33 

00 

00 

00 

]  00 

00 

33 

66  =  (11  -  12)  /  2 


Matrix  entries  show  only  subscripts. 

TABLE  5.  ELASTOPIEZODI ELECTRIC  MATRICES  FOR  POINT  GROUP  3m: 


THE 

[d]  , 

[g] , 

AND  [b] 

SETS. 

11 

12 

13 

14 

00 

00  ] 

00 

-22 

31 

] 

SE 

]  d ' 

12 

11 

13 

-14 

00 

00  ] 

00 

22 

31 

— 

-] - 

] 

d 

] (eps)T 

13 

13 

33 

00 

00 

00  ] 

00 

00 

33 

14 

-14 

00 

44 

00 

J 

00  ] 

00 

15 

00 

] 

sD 

]  g' 

00 

00 

00 

00 

44 

14*2  ] 

15 

00 

00 

— 

-3 - 

] 

g 

] (bet) T 

00 

00 

00 

00 

14*2  66  ] 

-22*2 

00 

00 

00 

00 

00 

00 

15 

J 

-22*2  ] 

11 

00 

00 

] 

SP 

3  b’ 

-22 

22 

00 

15 

00 

00  ] 

00 

11 

00 

— 

■3 - 

] 

b 

] ( zet) T 

31 

31 

33 

00 

00 

00  ] 

00 

00 

33 

66  =  (11  -  12)  *  2 


Matrix  elements  show  only  subscripts. 
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C3 3  =  (Sll  +  sl2 )  /  A 

C4 4  =  (Sll  -  Sl2)  /  B 

c66  =  (ell  -  cl2 )  /  2  =  s44  /  (2  *  B) 

K  =  c3 3  *  (ell  +  Cl2 )  -  2  *  cl3  *  Cl3 

L  =  c4 4  *  (Cll  -  Cl2 )  -  2  *  Cl4  *  Cl4 

2  *  sll  =  c33  /  K  +  c44  /  L 

2  *  Sl2  =  c33  /  K  -  c4 4  /  L 

sl3  =  -  Cl3  /  K  ;  sl4  =  -  cl4  /  L 

s33  =  (cll  +  cl2 )  /  K 

S44  =  (Cll  -  Cl2 )  /  L 

s66  =  (sll  -  sl2 )  =  2  *  c44  /  L 

det  (4x4)  [s]  =  A  *  B 

det  (4x4)  [c]  =  K  *  L 

A*K  =  B*L  =  A*B*K*L=1 


(80) 

(81) 

(82) 

(83) 

(84) 

(85) 

(86) 

(87a),  (87b)  (87) 

(88) 

(89) 

(90) 

(91) 

(92) 

(93) 


Formulas  (75)  to  (93)  hold  for  each  set  of  constant  electrical 


conditions: 

either 

E, 

D, 

or  P  constant. 

[CD]  ■ 

-  [cE] 

= 

[del 

cDE  ] 

= 

[e] •  [h]  = 

[h] •  [e] 

(23) 

del 

cDEll 

= 

+ 

e22 

h22 

+ 

e3 1  h3 1 

(94) 

del 

CDE12 

= 

- 

e22 

h22 

+ 

e3 1  h3 1 

(95) 

del 

CDE13 

= 

+ 

e31 

h33 

= 

+  h3 1  e3  3 

(96) 

del 

CDE14 

= 

- 

e22 

hl5 

= 

-  h22  el5 

(97) 

del 

CDE3  3 

= 

+ 

e3  3 

h3  3 

(98) 

del 

cDE4  4 

= 

+ 

el5 

hl5 

(99) 

del 

CDE66 

= 

+ 

e22 

h22 

(100) 

[CP]  ■ 

-  [cD] 

= 

[del 

cPD] 

= 

[a]  '  [h]  * 

( eps)  o 

= 

[h] •  [a]  * 

( eps)  o 

(24) 

del 

cPDll 

= 

( 

+  a2 2  h2 2 

+  a3 1  h3 1  ) 

*  (eps)o 

(101) 

del 

cPDl  2 

= 

( 

-  a22  h22 

+  a  3 1  h3 1  ) 

*  (eps)o 

(102) 
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del  CPD13 


=  (  +  a3 1  h3 3  )  *  (eps)o 

=  (  +  h31  a33  )  *  (eps)o  (103) 

del  CPD14  =  (  -  a22  hl5  )  *  (eps)o 

=  (  -  h22  al5  )  *  (eps)o  (104) 

del  CPD33  =  (  +  a3 3  h33  )  *  (eps)o  (105) 

del  cPD44  =  (  +  al5  hl5  )  *  (eps)o  (106) 

del  CPD66  =  (  +  a22  h22  )  *  (eps)o  (107) 

[cP]  -  [cE]  =  [del  cPE ]  =  [e]'  [a]  =  [a]'  [e]  (25) 

del  cDEll  =  +  e22  a22  +  e31  a31  (108) 

del  CPE12  =  -  e22  a22  +  e31  a31  (109) 

del  CPE13  =  +  e3 1  a33  =  +  a31  e33  (110) 

del  CPE14  =  -  e22  al5  =  -  a22  el5  (111) 

del  CPE33  =  +  e3 3  a33  (112) 

del  CPE44  =  +  el5  al5  (113) 

del  cPE66  =  +  e22  a22  (114) 

From  the  del  cl3  entries  we  have  the  ratios 

e3 1  /  e3 3  =  h31  /  h33  =  a31  /  a33.  (115) 

From  the  del  cl4  entries  we  have  the  further  ratios 

el5  /  e22  =  hl5  /  h22  =  al5  /  a22.  (116) 

[sE]  -  [sD]  =  [del  sED]  =  [d] '  [g]  =  [g] '  [d]  (26) 

del  sEDll  =  +  d2 2  g22  +  d31  g31  (117) 

del  SED12  =  -  d2 2  g2 2  +  d31  g31  (118) 

del  sEDl 3  =  +  d3 1  g33  =  +  g31  d33  (119) 

del  SED14  =  -  d2 2  gl5  =  -  g22  dl5  (120) 

del  SED33  =  +  d33  g33  (121) 

del  SED44  =  +  dl5  gl5  (122) 

del  SED66  =  +  d2 2  g22  *  4  (123) 
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[sD]  -  [sP]  =  [g]'  [b]  *  (eps)o 

=  [b]'  [g]  *  (eps)o  (27) 

del  sDPll  =  (  +  g22  b22  +  g31  b31  )  *  (eps)o  (124) 

del  SDP12  =  (  -  g22  b22  +  g31  b31  )  *  (eps)o  (125) 

del  sDP13  =  (  +  g3 1  b33  )  *  (eps)o 

=  (  +  b31  g33  )  *  (eps)o  (126) 

del  SDP14  =  (  -  g22  bl5  )  *  (eps)o 

=  (  -  b22  gl5  )  *  (eps)o  (127) 

del  sDP33  =  (  +  g3 3  b33  )  *  (eps)o  (128) 

del  sDP44  =  (  +  gl5  bl5  )  *  (eps)o  (129) 

del  sDP66  =  (  +  g22  b22  )  *  4  *  (eps)o  (130) 

[sE]  -  [sP]  =  [del  sEP]  =  [ b ] •  [d]  =  [d] '  [b]  (28) 

del  sEPll  =  +  d22  b22  +  d31  b31  (131) 

del  SEP12  =  -  d22  b22  +  d31  b31  (132) 

del  SEP13  =  +  d3 1  b33  -  +  b31  d33  (133) 

del  SEP14  =  -  d22  bl5  =  -  b22  dl5  (134) 

del  SEP33  =  +  d3 3  b33  (135) 

del  SEP44  =  +  dl5  bl5  (136) 

del  SEP66  =  +  d22  b22  *  4  (137) 

From  the  del  sl3  entries  we  have  the  ratios 

d31  /  d3  3  =  g3 1  /  g3  3  =  b31  /  b33  .  (138) 

From  the  del  sl4  entries  we  have  the  further  ratios 

dl5  /  d22  =  gl5  /  g22  =  bl5  /  b22.  (139) 

Piezoelectric: 

[d]  =  fe]  [sE]  (33) 

dl5  =  +  el5  sE4 4  -  e2 2  sE14  *  2  (140) 

d2 2  =  +  e22  ( sEl 1  -  sE12)  -  el5  sE14  (141) 
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d31  =  +  e3 1  (sEll  +  sE12 )  +  e33  sE13  (142) 

d3 3  =  +  e3 3  sE33  +  el3  sE13  *  2  (143) 

[h]  =  [ (bet) S]  [e]  (34) 

hl5  =  (bet) SI 1  el5  (144) 

h22  =  (bet)Sll  e22  (145) 

h31  =  (bet) S33  e31  (146) 

h33  =  (bet) S33  e33  (147) 

[g]  =  [ (bet) T]  [d]  (35) 

gl5  =  (bet) Til  dl5  (148) 

g22  =  (bet) Til  d22  (149) 

g31  =  (bet) T33  d31  (150) 

g3 3  =  (bet ) T3 3  d33  (151) 

[a]  =  [ ( zet ) S ]  [e]  (36) 

al5  =  (zet)Sll  el5  (152) 

a22  =  (zet)Sll  e22  (153) 

a3 1  =  ( zet) S33  e31  (154) 

a33  =  ( zet) S33  e31  (155) 

[b]  =  [ ( zet) T]  [d]  (37) 

bl5  =  ( zet) Til  dl5  (156) 

b22  =  ( zet) Til  d22  (157) 

b3 1  =  ( zet ) T3 3  d3 1  (158) 

b3 3  =  (zet ) T3 3  d33  (159) 

Dielectric : 

[ (bet) Y ]  =  ( (eps) Y] (-1)  (21) 

(bet)Yll  =  1  /  (eps)Yll  (160) 

(bet) Y33  =  1  /  (eps) Y33  (161) 
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(162) 


[ ( zetr ) Y ]  =  [ (betr) Y]  [I  -  (betr)Y](_1) 

(zet)Yll  =  1  /  ((eps)Yll  -  (eps)o)  (163) 

(zet)Y33  =  1  /  ( (eps) Y33  -  (eps)o)  (164) 

[(eps)T  -  (eps) S ]  =  [del  (eps)]  =  [e]  [d] '  = 

((chi)T  -  (chi)S]  =  [del  (chi)]  =  [d]  [e] '  (30) 

del  (eps) 11  =  del  (chi) 11  =  +  el5  dl5  +  e22  d22  *  2  (165) 

del  (eps) 33  =  del  (chi)33  =  +  e33  d33  +  e31  d31  *  2  (166) 

[(bet)S  -  (bet) T]  =  [h]  [g] •  =  [g]  [h] •  (31) 

del  (bet) 11  =  +  hl5  gl5  +  h22  g22  *  2  (167) 

del  (bet) 33  =  +  h33  g33  +  h31  g31  *  2  (168) 

[ ( zet) S  -  (zet)T]  =  [del  (zet)]  =  [a]  [b] '  =  [b]  [a]'  (169) 

del  (zet) 11  =  +  al5  bl5  +  a22  b22  *  2  (170) 

del  ( zet) 33  =  +  a33  b33  +  a31  b31  *  2  (171) 


16 


INPUT  VALUES  FOR  LI  TA  03 

The  values  measured  by  Smith  and  Welsh  (Ref.l)  using  the  pulse- 
echo  (transit-time)  technique  are  as  follows: 


TABLE  6 

cEll 

.  ISAGRIC  ELASTIC 

STIFFNESSES. 

cE12 

cE13 

cE14  cE3  3 

cE44 

CE66 

229.8 

Units: 

44 . 0 

10'  9) 

81.2 

N/m2 . 

-10.4  279.8 

96.8 

92.9 

TABLE  7 

PIEZOELECTRIC  STRESS  CONSTANTS. 

el5 

e22 

e3 1 

e3  3 

2.72 

1.67 

-0.38 

1.09 

Units: 

C/m2  . 

TABLE  8.  DIELECTRIC  PERMITTIVITIES  AT  CONSTANT  STRAIN. 
(eps)Sll  (eps)S33 

377.  379. 

Units:  10(-12)  F/m. 


OUTPUT  VALUES  FOR  LI  TA  03 

The  input  values  from  Tables  6,  7,  and  8  were  used  to  compute  the 
remaining  elastic,  piezoelectric,  and  dielectric  quantities  for 
lithium  tantalate  in  the  manner  discussed  in  prior  sections  of  this 
report.  The  results  are  given  in  Tables  9  to  16. 
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TABLE  9.  ELASTIC  STIFFNESSES. 


CE 

cD 

CP 

del  cDE 

del  cPE 

del  cPD 

11 

229.8 

237 . 6 

237.8 

7.78 

7.976 

0.187 

12 

44.0 

37.0 

36.8 

-7 . 02 

-7.19 

-0.169 

13 

81.2 

80.1 

80.1 

-1.09 

-1. 12 

-0.0261 

14 

-10.4 

-22.4 

-22.7 

-12.0 

-12.3 

-0.291 

33 

279.8 

282.9 

283 . 0 

3 . 13 

3.21 

0.0750 

44 

96.8 

116.4 

116.9 

19.6 

20.1 

0.472 

66 

92 . 9 

100. 3 

100.5 

7.40 

7.58 

0.178 

Units:  10 (  9)  N/m2. 


TABLE  10.  ELASTIC  COMPLIANCES. 


sE 

sD 

del  sED 

del  sEP 

del  sDP 

11 

4.93 

4 . 79 

4.78 

0.143 

0. 146 

0.00283 

12 

-0.518 

-0.424 

-0.422 

-0.0945 

-0.0962 

-0.00169 

13 

-1.28 

-1.24 

-1.23 

0.0450 

-0.0461 

-0.00105 

14 

0.585 

1 . 00 

1.01 

-0.419 

-0.427 

-0.00798 

33 

4 . 32 

4 .23 

4 .23 

0.0832 

0.0852 

0.00195 

44 

10.46 

8.98 

8.95 

1.48 

1.51 

0.0282 

66 

10.90 

10.42 

10.41 

0.475 

0.484 

0.00905 

Units : 

10 (-12) 

m2/N , 

TABLE  11. 

PIEZOELECTRIC 

[e] 

,  [h],  AND  [a] 

VALUES . 

e 

h 

a 

15 

2.72 

7.21 

7 . 39 

22 

1.67 

4.43 

4 . 54 

31 

-0 . 38 

-1.00 

-1.03 

33 

1.09 

2.88 

2.94 

Units:  e: 

C  /  m2 ;  h  and 

a : 

1q(  9)  V/m. 
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TABLE  12. 

PIEZOELECTRIC  [ d ] ,  [g],  AND 

[b] 

VALUES. 

d 

g 

b 

15 

26.5 

55.9 

56.9 

22 

7 . 51 

15.8 

16.1 

31 

-3 . 07 

-7.93 

-8.11 

33 

5.68 

14.7 

15.0 

Units:  d: 

10('12)  m 

/  V;  g  and  b:  10 

(- 

3)  m2/C. 

TABLE  13. 

DIELECTRIC 

(eps)  VALUES. 

(eps)  S 

(eps)T 

del  (eps)TS 

11 

377  . 

474  . 

97 . 1 

33 

247 . 0 

254. 6 

8.52 

Units:  lot 

-12)  F/m. 

del  (eps)TS 

=  del  (chi 

)  TS 

TABLE  14. 

DIELECTRIC 

(chi)  VALUES. 

II  l 

li  i 

il  l 

II  l 

ll  l 

il  l 

il  i 

il  l 

il  i 

n  l 

il  i 

(chi)S 

(chi)T 

i  li 

i  n 

i  il 

del  (chi ) TS 

11 

368  . 

465.7 

97.1 

33 

370. 

479  . 

8.52 

Units:  10 ( 

_12)  F/m. 

del  (chi) TS 

=  del  (eps)TS 

TABLE  15. 

DIELECTRIC 

(bet)  VALUES. 

(bet) S 

( bet )  T 

del  (bet) TS 

11 

2 . 65 

2.11 

-0.543 

33 

2 . 64 

2 . 58 

-0.0580 

Units:  10^ 

m/F. 
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TABLE  16.  DIELECTRIC  (zet)  VALUES. 


( zet) S 

( zet) T 

del  (zet)TS 

11 

2.72 

2  .  15 

-0.567 

33 

2.70 

2 . 64 

-0.0608 

Units: 

10(  9)  m/F. 

CONCLUSIONS 

This  report  provides  formulas  interrelating  the  coefficients  that 
appear  in  the  several  alternative  sets  of  constitutive  equations 
involving  the  elastic,  piezoelectric,  and  dielectric  properties  of 
crystals.  These  are  then  specialized  for  crystals  of  class  3m; 
using  measured  values  reported  for  lithium  tantalate,  numerical 
values  of  the  elements  of  the  polarization  matrices  are  calculated 
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483  Director 

US  Army  Material  Systems  Analysis  Actv 
ATTN:  DRXSY-MP 

001  Aberdeen  Proving  Ground,  MO  21005 

563  Commander,  AMC 
ATTN:  AMCDE-SC 
5001  Eisenhower  Ave . 

001  Alexandria,  VA  22333-0001 

609  Commander,  LABCOM 

ATTN:  AMSLC-CG,  CD,  CS  (In  turn) 

2800  Powder  Mi  1 1  Road 
001  Adel  phi ,  Md  20783-1145 

612  Commander,  LABCOM 
ATTN:  AMSLC-CT 
2800  Powder  Mill  Road 
C01  Adel  phi ,  MO  20783-1145 

680  Commander, 

US  Army  Laboratory  Command 
Fort  Monmouth,  NJ  07703-5000 

1  -  SLCET-00 

2  -  SLCET-DT  (M.  Howard) 

1  -  SLCET-DB 

35  -  Originating  Office 

681  Commander,  CECOM 

RSD  Technical  Library 

Fort  Monmouth,  NJ  07703-5000 

1  -  ASQfJC-ELC- 1 -T  (Tech  Library) 

3  -  ASONC-ELC- I -T  (STINFO) 

705  Advisory  Group  on  Electron  Devices 
201  Varick  Street,  9th  Floor 
002  New  York ,  NY  10014-4877 
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SUPPLEMENTAL  CONTRACT  DISTRIBUTION  LIST 
(ELECTIVE) 


205 

Di rector 

603 

Cdr,  Atmospheric  Sciences  Lab 

Naval  Research  Laboratory 

LAB COM 

ATTN:  COOE  2627 

ATTN:  SLCAS-SY-S 

001 

Washington,  OC  20375-5000 

001 

White  Sands  Missile  Range,  NM  88002 

221 

Cdr,  PM  JTFUSION 

607 

Cdr,  Harry  Diamond  Laboratories 

ATTN:  JTF 

ATTN:  SLCHD-CO,  TD  (In  turn) 

1500  Planning  Research  Drive 

2800  Powder  Mill  Road 

001 

McLean,  VA  22102 

001 

Adel  phi ,  MD  20783-1145 

301  Rome  Air  Development  Center 

ATTN;  Documents  Library  (TILO) 
001  Griffiss  AFB,  NY  13441 


437  Deputy  for  Science  S  Technology 
Office,  Asst  Sec  Army  (R&O) 

001  Washington,  DC  20310 

438  HQOA  (DAMA-ARZ-D/Or.  F.D.  Verderame) 

001  Washington,  DC  20310 

S20  Dir,  Electronic  Warfare/Reconnaissance 
Surveillance  and  Target  Acquisition  Ctr 
ATTN:  AMSEL-EW-D 

001  Fort  Monmouth,  NJ  07703-5000 

523  Dir,  Reconnaissance  Surveillance  and 
Target  Acquisition  Systems  Directorate 
ATTN:  AMSEL-EW-DR 

001  Fort  Monmouth,  NJ  07703-5000 

524  Cdr,  Marine  Corps  Liaison  Office 

ATTN:  AMSEL-LN-MC 

001  Fort  Monmouth,  NJ  07703-5000 

564  Dir,  US  Army  Signals  Warfare  Ctr 
ATTN:  AMSEL-SW-OS 
Vint  Hill  Farms  Station 
001  Warrenton,  VA  22186-5100 

602  Dir,  Night  Vision  £  Electro-Optics  Ctr 
CECOM 

ATTN:  AMSEL-NV-D 

001  Fort  Bel  voir ,  VA  22060-5677 


